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SUMMARY 


This paper presents a theoretical and experimental investigation 
of the effect of wind-tunnel walls on the air forces on an oscillating 
wing in two-dimensional subsonic compressible flow. A method of solving 
an integral equation which relates the downwash on a wing to the unknown 
loading is given, and some comparisons are made between the theoretical 
results and the experimental results. A resonance condition, which was 
predicted by theory in a previous report (NACA Rep. II50), is shown 
experimentally to exist. In addition, application of the analysis is 
made to a number of examples in order to illustrate the influence of 
walls due to variations in frequency of oscillation, Mach number, and 
ratio of tunnel height to wing semichord. 


INTRODUCTION 


In the evaluation of results obtained by measiarement of the forces 
on a wing in a wind tunnel, the question of the effect of the tunnel 
walls arises. In the case of steady flow the problem has been extensively 
investigated and, in general, relatively simple factors have been deter- 
mined which can be used to modify measurements of the forces on a wing 
in a tunnel to correspond to free-air conditions. However, the corre- 
sponding problem of the effect of walls on an oscillating airfoil has 
received relatively little attention, particularly in the case of com- 
pressible flow. The present paper concerns the wall effects in the 
oscillating case and treats the problem in two-dimensional subsonic 
compressible flow. 

Hi incon 5 )ressible flow, theoretical treatments of wall effects on 
oscillating wings have been made by several investigators and reported 

^Supersedes and extends recently declassified NACA RM L52I17a, 1955. 
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in references 1, 2, and 5* These investigators have shown generally 
that the tunnel-wall effects are a maximum for some small values of the 
reduced frequency and that the wall effects become negligible as the 
reduced frequency is increased. Extension of the theoretical treatment 
of the problem to include the effects of con^jressibility of the fluid 
has been reported in reference k. In this reference, it is shown that, 
in addition to the large effect noted at low values of the reduced fre- 
quency, imder certain conditions, large effects of the walls may be 
encovmtered at hi^er values of the reduced frequency. These effects 
are due to an acoustic resonant phenomenon which occurs when a disturbance 
from the oscillating wing is reflected from the txmnel wall back to the 
wing with such a phase relationship that it reinforces a succeeding 
disturbance. 

In reference 4 the problem was expressed as an integral equation 
which relates the known downwash distribution over the airfoil to the 
Tinknown lift distribution. One purpose of the present paper is to dis- 
cuss further the integral equation and to demonstrate a method of solving 
it. A second pxirpose is to present some results showing wall effects 
calculated by this procedvire and, in some cases, to compare the calcu- 
lated results with experimental results. This phase of the investigation 
is given in three parts: (l) a comparison between analytically and exper- 

imentally determined values for the lift and moment on a wing oscillating 
in pitch at several subsonic Wh.ch numbers j (2) an analytical study of the 
effects of a variation in Mach number for a constant ratio of tunnel height 
to wing semichordj suid ( 5 ) an analytical study of the effects of a vari- 
ation in the ratio of tunnel hei^t to wing semichord. Portions of this 
material have been reported previously in reference 5 QJ^d are included in 
the present paper in order to provide a more extensive and unified 
presentation. 

As a check, the integral equation for the downwash on a wing oscil- 
lating between walls in a compressible medium is reduced to the zero- 
frequency condition and is shown in an appendix. The resulting expression 
is in agreement with steady-state results. 

The calcvilation procedxare and the results contained in this paper are 
of significance for such problems as the experimental measurement of the 
forces on em oscillating airfoil, the determination of wing-flutter char- 
acteristics in wind txmnels, and also in certain possible types of flutter 
of airfoils in cascade. 


SYMBOLS 


a velocity of soimd, ft/sec 

Ajj coefficients in series expression for lift distri- 

bution (eq. ( 16 )) , where n = 0 , 1 , 2 , . . . 
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b 

h 

1 

H 

Ho(2) ,Hi(2) 

k 

K(M,z) + K(M,z,H) 

L(xo),L(9o) 

lu 

Hi 

Hi 

M 


■D = 

2rt3 


ving semi chord, ft 

displacement of wing in vertical translation, ft 
height of tunnel, ft 

height of tunnel referred to wing semichord 
Hankel functions of the second kind 

reduced- frequency parameter, boj/U 

kernel of integral equation 

lift distribution, Ib/ft/unit span 

aerodynamic lift force per unit span due to pitch 

eterodynamic lift force per unit span due to translation 

aerodynamic moment per unit span due to pitch 

aerodynamic moment per unit span due to translation 

Mach number, U/a \ 


Rn = \/(x - Xo)2 + 
U 

w(x) 

% 

x,Xo,y,| 

z = k(x - xo) 

a 

3 = )ll - 


3^(nH)^, where n = 1,2,5, • • • 

stream velocity in chordwise direction, ft/sec 

vertical induced velocity (perpendicular to chord), 
ft/sec 

axis of rotation measured from midchord, positive 
rearward, based on semichord 

Cartesian coordinates 


angular displacement of wing in pitch, radians 
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e = | x - Xq I 

pH 

9 = -cos"^ 

0Q = -COS“^Xq 

fluid density, slugs/cu ft 

phase angle hetween lift force and position of 
pitching wing, deg 

phase angle hetween lift force and position of 
translating wing, deg 

phase angle hetween moment and position of pitching 
wing, deg 

phase angle hetween moment and position of trans- 
lating wing, deg 

circular frequency of oscillation, radians/sec 

(Ores circular frequency at resonance, radians/sec 

Ap pressure difference hetween upper and lower surface, 

Ih/sq ft 

Primed qviantities refer to a wing in free air. 


K 

K 

Hi 


ANALYTICAL INVESTIGATION 


This section is concerned with the development of a method for 
solving the integral equation, originally derived in reference 4 , which 
relates the downwash to the loading on an oscillating wing. The basic 
integral eq\iation and its kernel is given hy equations (l) and (2). 
Reduction of the kernel is made in equations (5) to (lO). Alternative 
series expressions for the kernel which are suitable for numerical 
con5)utation are given hy equations (U) to (15) . The loading on the 
wing is given hy equation (16), the downwash expression hy equations (18) 
and (19) , and finally the lift and moment expressions hy equations ( 20 ). 
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The Integral Equation and Its Kernel Function 

The integral equation .- The integral equation of reference 4 for 
the vertical velocity or downwash of an oscillating airfoil between 
plane walls may he written as 


w(x) - ^^^2 I'(^o) 2) + dxo 


( 1 ) 


where w(x) is the known vertical velocity (or known motion of the wing) 
and L(xo) is the unknown lift distribution or the local strength of a 
distribution of oscillating pressure doublets. The functions within the 
brackets con5)rise the kernel function of the integral equation and appear 
formally as 


k(M,z) 


^ lim e-ik(x-xo) 



(2a) 


K(M,z,H) = 11m f 

V-^n 


— P 

2 Z (-l)" & Ho<2) 

n=l 


^ \/|2 + p2(y 


nH)2 


d| (2b) 


The first function K(M,z) corresponds to the kernel for the free-air 
condition as given by Possio (ref. 6). The second function K(m,z,H), 
containing the infinite summation, is the additional part of the kernel 
arising from the effect of the walls. Physically, a kernel function 
represents the contribution to the vertical velocity at a field point 
due to a pulsating pressure doublet of unit strength located at any other 
point in the field. For the particular case represented by equations (2), 
the kernel function gives the vertical velocity in the plane of a wing 
located in the center of the tunnel. The expression K(M,z) gives the 
downwash of a doublet in the plane of the wing, whereas the eiqoression 
K(M,z,H) gives the downwash due to the system of images which mathemat- 
ically represents the walls. 

Reduction of the kernel function .- The integrals contained in the 
expressions for the kernel function in equations (2) are improper because 
they have an infinite limit and also because, at certain points, the 
integrands become singular. This section is concerned with the reduction 
of these integrals to a form more amenable to computation. 

By making use of the fact that the Hankel functions in equations (2) 
satisfy the identity 
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The integrals of equation (l^•) that contain partial derivatives of 
Hankel functions can he integrated twice by parts to obtain 
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The last integral of equation ( 5 ) may be written in two parts as 



The first integral on the right-hand side of equation (6) will be left 
temporarily in integral form and will be treated in the following section, 
(See evaluation of S5 following eq. (I 5 ).) 

The second integral on the right-hand side of equation (6) has not 
been integrated in closed form) however, in wind-tunnel problems it can 
be handled conveniently by approximate methods. (An alternative means 
of treating this integral, which avoids the approximation, but is somewhat 
more tedious, will be indicated in the section following eq. 14(c).) 

A practical assumption which is often made in the analysis of the effect 
of wind-tunnel walls is that the tvinnel height is considered large compared 
with the wing semichord. With this assumption the argimient of the Hankel 
function in equation (6) can be written as (in the limit as y — »0) 


provided that — ^ « 1. 

PnH, 

This approximation implies that the airfoil images, and, in partic- 
ular, the closest image n = 1, are a sufficient distance from the airfoil 

so that the actual distance \/|2 + 32(ng) ^ may be replaced by the vertical 
distance pnH of the image above the airfoil. Of course, this approxima- 
tion does not hold for Mach nvunbers close to or equal to \mity. With this 
approximation, the second integral of eqiiation (6) can be expressed as 
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and these equations may he used to express equation (4) as 


w(x) 



+ 


K(M,z,H) 


dxo 


(7) 


( 8 ) 


where 


K(m,z) = ^ e"^^ 
4p 


I2 

eP^ 


-Ho(2)(^Ro) + Hi(2)(^Ro) 

|X - Xo| 


ip" 


1 +P ^ / M ^iu 


M 


— logg ^ ^ ^ ^ Ho^ 2) (M(u|)(iu 

rtP 


(9) 


and 


K(M,z,H) = ^ 


-iz 


iz 


-eP2 21 (-I)"" Ho(2)(nRn) + 


n=l 


(eS2 - 1 ) ^ (- 1 )" 


00 nco ^ i.kfe 

ik^ (-1)^ / e*P^ 
n=l ^0 


Hb 


(2) 


M/|2 + p2(njj)2 


d.^ + 


■Wl 


n=l 




( 10 ) 
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in which use has been made of 







Ro = |x - Xq 


R, 


n 


/(x - Xo)^ + 32(nH)2 


Equation ( 8 ), together with the definition of equations (9) and (lO), 
permits the determination of the effect of tunnel walls on a lift distri- 
bution L( 9 o) for a given downwash distribution w(x) . The integral 
equation for the case of no tunnel walls checks the results of Possio 
(ref. 6 ) . For the case with walls and for the limiting steady-flow case 
of zero frequency, it is possible to obtain a mathematical check with 
some existing results; this check is shown in the appendix. 


Alternative series expressions for kernel .- Although the fonn of the 
kernel K(M,z,H), given by equation (lO), could be used for calculation, 
alternative series which are more hi^ly convergent may be used and are 
given in this section. 

The kernel K(M, z,H) is the s\jm of four infinite series which can 
be written as 


K(M,z,H) = + CgSg + C5S5 + C4S4) ( 11 ) 

where the S^'s denote the indicated infinite summations of equation (lO) 
and the Cn*s the respective multipliers. 

\ 

Series Sp and S2 of equation (ll) may be put in a more rapidly 

convergent form according to Infeld, Smith, and Chien (ref. j) . When the 
variables p and e are introduced, where 


and 


^ 2 np 



the series S^ and S2 can be written as 
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Si =2l (-I)"" 

n=l 


n=l 


_1_ 

2rt 


(-1)" Ho^^^(2np l/e2 + n2 ) 

Tie }}( 2n+l) 2-i|.p2 


^ + 2i > ^ 


|_ \/l - 4 p 2 n=l ^(2n + l)2 - l^p2 


e-«€>/(gn-l)£-4p2 _ 

\/(2n - 1)2 - 4 p 2 


and 


S2 (-1)"" 

n=l \ P / 


(-1)“ Ho^2)(2jtnp) 


n=l 


2it 


( 12 ) 


00 

-It + 21 (y + loffo / 

N 

1 1 ~| 

n=l 

^/(2n - 1)^ - 2n - 1 


( 15 ) 


where Euler's constant 7 = 0 . 577215 * 

Series S5 may he evaluated hy utilizing the expression for Si 
(eq. (12)) and integrating the resulting expression to obtain 


®5 ", 


r e‘p2^ 

00 


21 (-1)'' Hq(2) 

^ \/|2 + p2(nH)2 

'0 

n=l 



d| 


(l 4 a) 
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= - A log^ 
rtk 


1 + B + 

M 


21 B ^ 1 M_ 

“ - (-/[(- . - V] . 

(I 4 c) 

It is of interest to note that series Sj may he employed in, an 

alternative means of integrating equation (7). For application to wind- 
t\mnel problems, where the ratio of timnel hei^t to wing semichord is 
small, or in application to cascade problems, the approximation employed 
in integrating equation (7) becomes less valid. It is possible to avoid 
the use of the approximation by writing the integral of equation (7) in 
a form which is identical to that of equations (l 4 a) and (l 4 b) with the 
exception of the upper limit. The integral containing the Hankel function 
can be evalmted by employing the tables of Schwarz (ref. 8 ) . The second 
integral, containing only an exponential term, can be integrated in 
closed form, as was done to obtain equation (ikc). 

Series Sij. may be evaluated in a direct manner by employing tables 
of the Hankel function and by using for large values of the argument the 
approximation 

Hi(2)(^Rn) « 

With the aid of series S^, S2> S5, and S4, the kernel K(M,z,H) may 

be evaluated. 
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Method of Solution 

A method of using equation (8) to determine the aerodynamic forces 
on a wing oscillating in the presence of plane walls is now discussed. 

The method under consideration is one of collocation similar to that 
used hy Possio (ref. 6) and Frazer (ref. 9) for the case of no walls. 

The approach involves the assumption of an' appropriate series expression 
for the lift distribution, substitution of this series in the integral 
equation for the downwash, and calculation of the downwash at arbitrarily, 
selected points on the chord (control points). Thus equation (8) is 
reduced to a set of simultaneous equations, the unknowns of which are 
the coefficients of the assumed expression for the loading. 

Expression for the loading . - The expression which is assumed for the 
lift distribution is a trigonometric series expansion which satisfies the 
Kutta condition at the trailing edge and which has the proper type of 
singularity at the leading edge. This expression is 


00 " 

HigI = Aq cot -^ + ^ An sin n9o = L(0o) (16) 

pU^ ^ n=l 

where Xq = -cos 9 q and the An's are unknown coefficients to be deter- 
mined in accordance with the downwash w(x), which is known from the 
motion of the wing. It is desirable to rewrite equation (8) in terms of 
the variable 9 q as follows: 


w(x) = Uk 



L(9o) K(M,z) sin 9^ d9o + 



L(9o) K(M,z,H) sin 9 q d9o 


(17) 


The first integral on the right-hand side of equation (17) is the 
integral expression first derived by Possio (ref. 6) for the condition 
of no walls. Its solution has been treated by several investigators and 
will not be discussed here. It can be expressed entirely in terms of 
the unknown coefficients of equation (l6) . The second integral of 

equation (I7) may he evaluated hy the use of equations (12), (l5)^ 
and (15). 


Determination of the aerodynajuic forces . - The integrals of equa- 
tion (^IyT determined for a selected numher of control points and 
equated to the expression for the downvash. The expression relating the 
downwash to the motion of a wing translating (h) and pitching (a) about 
an axis located at x^ is 


w(x) = h + Ua + h(x - Xa)d 


(18) 
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or, with the assumption of harmonic motion. 



Equation (I 9 ) is used to calculate w(x) for values of x appropriate 
to each of the selected control points. A set of simultaneous equations 
can then he written, the number of which corresponds to the number of 
control points employed and (conveniently) to the number of terms 
retained in the series for L(9o)- The unknown coefficients may now be 
determined by solving these simultaneous equations. The total lift and 
moment about the midchord are given in terms of the coefficients ^n 
through the relations 

^ - |(jk) + I 

Effect of the number of control points considered .- An investigation 
was made of the number of terms of the series for the lift distribution 
(eq. ( 16 )) and thus of the number of control points required to obtain 
satisfactory accuracy. Calculations were performed for a particular case 
by increasing the number of control points and the number of terms of the 
loading series until the solutions were in reasonable agreement. For the 
case considered, three terms of the series for the lift and three control 
points at the quarter-, half-, aind three -quarter -chord positions gave 
satisfactory results. The consideration of two additional control points 
at the leading and trailing edges, together with two additional terms of 
the lift series, made no significant change in the results. For high 
values of the reduced-frequency parameter k, the use of additional 
control points might be necessary. 

The procedure just discussed involves consideration of a continuous 
distribution of pressure doublets over the chord. Calculations requiring 
much less computing can be made by considering the chordwise loading to 
be concentrated in a single doublet located at the quarter chord and 
satisfying the downwash at the three-quarter chord. In the case of the 
lift, this approach has been found to give fairly good agreement with 
the results of the more elaborate calculations except in the vicinity 
of the resonant frequency. 
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The Analytically Indicated Resonance Phenomenon 

Tvfo-dimensional tunnel .- By examination of eqiiations (12) and (15), 
it may be seen that the series Si and S 2 become infinite when 

= (2n - 1)2 

or where 

^ = rtp(2n - 1) (n = 1,2,5, • • •) (2l) 

cl 

At these critical values of the frequency parameter, the expression for 
the kernel K(M,z,H) (eq. (ll)) becomes infinite for all values of x. 
Physically, this condition represents a resonance in the tunnel involving 
a transverse oscillation of the moving air between the walls. 

The fundamental or smallest critical values of coH/a corresponding 
to n = 1 in equation (21) are shown plotted as functions of Nfeich 
number M in figvire 1. Equation (21) and figure 1 show that finite 
values of the critical frequency exist for the condition M = 0, U = 0, 
and a ^ 00 . These conditions correspond to a compressible fluid at zero 
velocity in the tunnel. As the Mach mamber is increased, the critical- 
frequency parameter decreases rapidly euid becomes zero at a Mach number 
of unity. 

As indicated by equation (l), the product of the lift and the kernel 
function must remain eqml to the vertical velocity over the wingj this 
velocity is defined by the motion of the wing and remains finite. The 
product of the lift and the kernel function can remain finite only if the 
lift approaches zero as the kernel becomes infinite. This condition in 
the tunnel is analogous to the well-known case of a simple undan^jed- 
spring-mass system for which, at the resonant frequency, theory predicts 
an infinite deflection of the mass occurring even with a forcing function 
of small amplitude. 

Circ;ilar tunnel .- A resonance can also be demonstrated for the 
infinite circular tunnel. The nature of the boundary-value problem, for 
this case, makes it possible to separate variables; therefore, the 
governing partial-differential equation can be reduced to Bessel’s 
equation. (See, for instance, ref. 10.) The resonant frequencies are 
then fovind as the roots of the equation 




(n ~ • . •) 
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or 


gJ) 

2a 


Pn3 


where Jjj represents the Bessel function of the first kind, D is 
the txinnel diameter, and Pj^ is the root of the equation 

Jn'(Pn) = 0 

Values for for the first several modes are Pj^ = 1.84, 5*05^ and 

4.17- Note that, for a circular tionnel having a diameter equal to the 
height of a plane tunnel, the fundamental frequency is 5*68/jt = l.l? 
higher than resonant frequency in the plane tunnel discussed in this 
paper . 


EXPERIMENTAL INVESTIGATION 
Wind Tunnel 


The experimental part of the investigation of the effect of tunnel 
walls on the forces acting on an oscillating airfoil was conducted in 
the langley 2- by 4-foot flutter research tiinnel. For these tests, a 
rectangular test section having dimensions of 2 feet by 5-8 feet was 
used. This tunnel is of the closed-throat, single-return type and 
en5)loys either air or Freon-12 as a testing medium at pressures from 
1 atmosphere down to about I/8 atmosphere. 

It has been shown previously that the resonant frequency varies 
directly as the si>eed of sound. Inasmuch as Freon- 12 has a speed of 
sound equal to about one-half that of air, the experiments to be dis- 
cussed were conducted in Freon-12 so that the resonant frequency could 
be surveyed within the frequency limitations of the equipment. 


Model and Oscillating Mschanism 

Figure 2 is a schematic drawing of the test section with the model 
eind oscillating mechanism installed. The model had a chord of 1 foot 
and MCA 65-OIO airfoil sections j it con^jletely spanned the 2-foot dimen- 
sion of the test section. The gaps between the model and the tiinnel 
wall were sealed by end plates which rotated with the model. The model, 
driven symmetrically from both ends, was oscillated in pitch about the 
midchord by a direct-drive eccentric-cam system powered by an induction 
motor with variable frequency supply. 
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Instrimentation 

The lift and moment on the wing were obtained by electrical integra- 
tion of the oul^juts of 12 model 49-TP NACA miniature electrical pressure 
gages. The pressure gages, which are described in considerable detail 
in reference 11, were located at the center of the span at 2.5> 5^ 10, 

15, 20, 50, 40, 50, 60, 70, 80, and 90 percent of the chord. Each gage 
was arranged to indicate the difference in pressure between orifices on 
the upper and lower surfaces. Electrical integration techniques used in 
these experiments are discvissed in reference 12. The so-called square- 
wave method of weighting was. used; that is, the pressure indicated by 
each gage was assumed to represent the pressure acting over a portion of 
the chord extending one-half the distance to the next gage both forward 
and rearward. For example, the fraction of the chord assigned to the 
first gage was 5*75 percent and to the sixth gage was 10 percent. Some 
of the implications of this method of integration will be discussed in 
a subsequent section. 

The angular displacement at the midspan position was indicated by 
resistance-wire strain gages attached to a torque rod running through 
the center of the hollow wing. One end of the torque rod was fixed to 
the center of the wing and the other end was fixed to the timnel wall. 

A schematic diagram of the instrumentation is shown in figure 5- 
The magnitude of the vector representing the fundamental conqjonent of 
lift or moment and angular displacement was indicated on an alternating- 
current vacuum-tube voltmeter attached to the output of a variable- 
frequency, narrow-pass-band filter. In essence, the filter performed 
the function of a Fourier analysis in that both random components and 
higher harmonics were removed from the signal. In order to measure the 
phase angle between lift or moment and the angular displacement, the 
output of the filter was fed into a pulse-shaping circuit designed to 
convert the sinusoidal signals into pulses corresponding in time to the 
"crossover" points of the original signal. The pulses were then fed 
into an electronic chronograph that accurately indicated the time 
interval' between the leading pulse which started the chronograph and 
the lagging pulse which stopped it. The ratio of this time Interval 
to the period of oscillation, when multiplied by 560°, yields the phase 
angle in degrees. The period and frequency of oscillation were deter- 
mined by starting and stopping the chronograph with the angular- 
displacement signal. In order to minimize the effects of small differ- 
ences in components between the two circuits, a "tare" switch was pro- 
vided which fed a single signal (the angular displacement) through both 
circuits. The resulting phase angle represented the phase shift intro- 
duced by the filters and pulse-shaping circuits. 
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Test Conditions 

The Mach number of the tests was varied from M = 0.^ to M = O.70 
and the Reynolds number was held constant at about 5 x 10° veirying the 
density. The frequency of oscillation was varied from 0 to 60 cycles per 
second, and the magnitude of angular displacement was about 1.2° except 
for some lift data at M = O.7I which was obtained at ein angular dis- 
placement of about 2.4°. 


DISCUSSION OF RESULTS 


The theory and calculation procedure and the experimental technique 
discussed previously for the determination of the forces acting on a 
wing oscillating between walls have been applied to a number of specific 
examples. The investigation has been divided into three parts: (l) A 

comparison is made of analytical and experimental results obtained for 
the lift and moment on a pitching wing for several subsonic Mach numbers, 
(2) theoretical results for the effects of a variation in Mach number at 
constant tunnel hei^t are given for a pitching wing and also for a wing 
undergoing vertical translation, and (5) theoretical results for the 
effects of a variation in the ratio of tunnel height to wing semichord 
are presented for particular values of Mach nvunber. 


Comparison of Theory and Experiment 

In figure 4 a con5)axison is made of analytical and experimental 
results for a wing oscillating in pitch about its midchord. Figures 4(a), 
4(b), 4(c), and 4(d) apply, respectively, to Mach numbers of 0.55# 0.5# 
0.6, and 0.7* The results apply to a ratio of tunnel height to wing 
semi chord H of 7*60. 

The plots on the left-hand side of each figure show the magnitudes 
of the forces and moments as a function of the frequency of oscillation 
whereas those on the right-hand side show the corresponding phase angles. 
The magnitudes are presented as ordinates in the form of ratios 
and 3h these ratios, the quantities and are, respec- 
tively, the lift force and the moment on a wing in a tunnel; and 

are the theoretical lift and the theoretical moment on a wing in 
free air. The effect of the tunnel walls appears, therefore, as a devia- 
tion from uniiy of the ratios and |Hx/Hx'| when Iq, and 

% are the theoretically derived forces and moments. When Iq, and % 
represent the experimental forces and moments, the deviation from unity 
may not be completely attributed to the effect of tunnel walls because 
such factors as airfoil thickness and viscosity may cause deviation from 
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the elementary theory. The abscissa in the figures is the ratio of the 
frequency of the pitching oscillation to a frequency calculated for the 
resonant condition. 

Excellent agreement between theory and experiment is obtained for 
the phase ajigles, in most cases, for both the lift and the moment. 
Quantitatively, however, the agreement between theory and experiment for 
the magnitudes of the forces is not as good, althou^ very similar trends 
are demonstrated; in most cases, a systematic difference appears. Some 
possible sources of the differences between theory and experiment are 
discussed in the following section. 

Examination of figure l4- reveals that the theory predicted the 
resonant frequency very well. In all cases, the minimum lift and 
moment was found to lie very close to the analytically indicated resonant 
frequency. Theoretically, the lift and moment reduce to zero at the 
resonant condition. Under actual conditions, such as finite tunnel 
length, transmission of energy through the walls, nonlinearities at 
higher aunplitudes, and turbulence in the flow that gives rise to damping, 
pure resonance is unobtainable. However, it may be seen by examining 
figure 4(d) that the lift and moment were reduced to 20 percent of the 
values away from resonance. 


Remarks on Some Differences Between Theory and Experiment 

In this section, some limitations and possible reasons for differ- 
ences between theory and experiment are listed and discussed. 

In the coD5)arison between theory and experiment shown in figure 4, 
an almost constant difference of 10 to I5 percent between the magnitiides 
is to be noted, whereas the phase angles are in good agreement. Of the 
several possible reasons for these differences between theory and exper- 
iment, perhaps the more important ones are airfoil thickness, Reynolds 
number, finite tunnel length, transmission of energy throu^ walls, 
dissipation of the pressure waves due to turbulence, and integration 
procedures. The effect of all these possibilities is not known for the 
oscillatory case. In the steady-state case, however, it is known that 
the effect of increasing thickness is to increase the slope of the lift 
curve . 

In the considerations of the analytical integration, a collocation 
scheme was used to solve the integral equation. In general, the three 
collocation points used were found to be satisfactory as pointed out 
previously. 

Twelve pressure cells were used for the experimental integration of 
the forces. A stepwise integration procedure was enployed; that is, the 
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pressure as recorded by the pressure cell was multiplied by an area of 
the wing over which it is assvimed that the cell will give an average 
pressure. This procedxire gives good results except possibly at the 
leading edge where the pressure variation is very great. Theoretically 

the pressxire approaches infinity at the leading edge ^as 

ejqperimentally it is found to be very large. As a matter of fact, if a 

theoretical distribution of pressure is assimed to be cot ^ = 

as shown in the following sketch ajid this curve is integrated in the 




Illustration of the differences in integration 


same manner as the e^erimental curve was integrated - that is, by calcu- 
lating the ordinate at the same value of the abscissa at which the 
pressure cells were located for the experiment - it is found that the 
area as determined by the approximate method is 8 percent less than the 

area as determined by integrating cot — in closed form. It is apparent 

that the neglected area (shaded) can be appreciable. In the actiial exper- 


iment, in which a fairly thick airfoil was used, the neglected area would 
probably be smaller but could perhaps contribute to the almost constant 
difference in magnitudes of lift and moment between the theoretical and 
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experimental results. The effect of this integration difference on the 
phase angles, which have heen shown to he in good agreement, would not 
he as pronoxinced. 


Effect of a Variation in Mach Nuniber at Constant Tunnel Height 

An emalytical investigation has heen made of the effects of a vari- 
ation in Mach number at constant tunnel height on the forces on an 
oscillating wing. Some of the results of the previous section are 
employed together with results of additional calculations . The magni- 
tudes and associated phase singles of both the lift force and the moment 
have heen determined for a pitching wing and also for a wing undergoing 
vertical translation. Calculations have heen made for a constant value 
of the height-semichord ratio H of J.SO eind for Ifech numbers of 0.3 > 
0.5, 0.7# and 0.8. Results of the calculations are shown in figures 5 
and 6 for the lift and in figures 7 and 8 for the moment. 


The magnitudes of the forces and moments are presented, as in the 
previous section, in the form of ratios: jLgyLg^'j, jMo,/l^'j, 

eind jMh/^h'l • phase angles related to these ratios are presented 


as a difference between a wing In a tunnel and a wing in free air. The 
magnitudes and phase angles are plotted against a frequency parameter 
o)H/a, where co is the circular frequency of oscillation of the wing, 

H is the height of the tunnel, and a is the velocity of sound. At a 
particular value of the frequency parameter, a progressively larger effect 
of the walls is indicated as the Mach number increases. At all Mach num- 
bers, the lift is reduced to zero at the resonant frequency. The dip in 
the curves against frequency ratio, which appears to he characteristic 
of the low Mach number cases, gradually disappears as the Mach number is 
increased. As in the case of the magnitude of the lift, the effect of 
the walls on the phase angle increases as the resonant frequency of the 
tunnel is approached and also as the Mach number is increased. 


The magnitude of the moment about the midchord is shown in figures 7 
and 8. The cvurves have the same shape as the lift-ratio curves and again 
decrease to zero at the resonant condition. The corresponding phase 
angles are shown in these figures. Note that in figures 5 to 8 only 
sligd^t differences appear between results for pitch and those for 
translation. 


Effect of Variation in Ratio of Tunnel Heigjht to Wing Semichord 

The effect on the lift-force ratio and on the associated phase angles 
of varying the ratio of tunnel hei^t to wing semichord is illiistrated in 
figure 9. The results presented in figures 7 and 8 have been based on the 
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consideration of a distribution of pressure doublets over the chord of 
the airfoil and on satisfying the downwash condition at three chordwlse 
stations. Results for figure 9 have been obtained by the simplified 
procedure of concentrating the loading in a single doublet at the 
qviarter chord and satisfying the downwash condition at the three-quarter 
chord. This approach gives fairly good agreement with the results of 
the more elaborate procedure except near the resonant frequency. 

Calculations have again been made for an airfoil oscillating in 
pitch about its midchord for values of the height-semichord ratio H 
of 7.60, 16, and 52 at M = O.5 and 0.8. In figure 9(a) the lift 
ratio |Iq,/Iq,'| is plotted against the reduced-frequency param- 
eter k = bo>/U. Plots for both ^fa.ch numbers are made to the same scale 
for ease of comparison. It is again apparent that the effect of 
reducing the Mach number is to reduce the effect of the tunnel walls 
and to raise the value of the critical frequency at which resonance 
can occur for a given tunnel. For example, for H = at M = 0.8, 

the critical value of k is O.5O, whereas for M = 0.3# the critical 
value of k is increased to I.3I. Also, as was to be expected, 
increasing the hei^t of the tunnel has a marked effect in reducing 
the influence of the tunnel walls for most of the frequency range. 
However, the critical frequency is reduced for increased tunnel height 
so that in the large tunnel the range of k below the fundamental 
resonance becomes smaller. This reduction in frequency would seem to 
be a disadvantage of the larger tunnels. However, the second branch of 
the curve for H = 32 at M = 0.3 shows that, for frequencies between 
the first and second resonant points, the effect of the walls on the 
magnitude of the lift is no greater than for frequencies below first 
resonance. Note also that the approach to resonance is quite abrupt; 
consequently, only a small range of frequencies very close to resonance 
would be critical and tests could then be conducted between the critical 
frequencies . 

In figure 9(^) "the phase angles, associated with the results of 
figure 9(a’) , for both the wing in a txmnel and for a wing in free air 
are shown as a fimction of reduced-frequency parameter k for values 
of hei^t-semichord ratio H of 7.60, 16, and 32 at M = 0.3 and 0.8. 

At M = 0.3 t the effect of walls on the phase angle is generally 
very small; at M = 0.8, the effect is small at low frequencies but 
increases greatly as the critical value of k is approached. As the 
resonant frequency is approached, results at both ^fe,ch numbers show that 
the phase angle increases in negative value and appears to approach -90°. 
As the resonant condition is exceeded, there is a sudden shift in phase 
angle. This change is similar to that found for the oscillation of a 
sinqple undan5>ed spring-mass system where an abrupt change in phase angle 
of 180° is found as the resonant frequency is exceeded. Because of the 
complexity of the kernel, which involves an infinite series of Hankel 
functions, the phase angle at resonance has not been determined. 
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In figure 10, the plot from figure 9(a) for M = 0.5 and H = 7-60 
is con 5 )ared with some results of Reissner (ref. 2) for the effect of 
walls on the lift-force ratio in incompressible flow. This curve for 
M = 0 significantly duplicates the rather large wall effect at low 
values of k which is noted for the M = 0.5 result. (At low values 
of k, the curves for the two different Mach numbers are almost coinci- 
dent probably because of the fact that the slightly lower height- 
semi chord ratio H = 6.28 used by Reissner counteracts the effect of 
decrease in Mach number.) For values of k greater than 0.5^ the 
cirrves for the two Mach niimbers separate; the lift-force ratio at 
M = 0 approaches unity and exhibits no effects of resonance because 
the resonant phenomenon arises only from the effects of con 5 )ressibility. 


CONCLUDING REMARKS 


This paper has dealt with the problem of an airfoil oscillating 
between plane walls in subsonic compressible flow. It constitutes a 
continiiation of the work initiated in NACA Rep. 1150 in that a method 
of solving the integral equation is presented and some experimental 
results are compared with theory. 

The comparison between theory and experiment for the phase angles 
between lift force or moment and position is shown to be very good, 
whereas the comparison between theory and experiment for the magnitudes 
of the lift and moment is not as good; however, the trends are all 
accxirately predicted. In all cases the resonant frequency was 
accurately predicted. The cause of the apparent discrepancy between 
the theoretical and experimental lift and moment may be attributed to 
several factors such as dissipation of the pressure waves due to turbu- 
lence of the air flow and transmission of the energy through the tunnel 
walls. In addition, the theoretical work was based on the concept of 
a very thin wing at Infinite Reynolds number, whereas the experiments 
were made with a 10-percent-thick wing at a Reynolds number of approxi- 
mately 5 X lo6. The effect of thickness and Reynolds number have not, 
as yet, been delineated for the oscillating case. 

As would be expected, it is shown theoretically that the larger 
the tunnel the less the effect of the walls. The critical frequency, 
however, is also reduced as the tunnel height is increased, but it is 
shown that tests may be made above the resonant frequency with no larger 
tunnel-wall effect than is found below the resonance. In addition, the 
range of influence of the resonant region is greatly reduced so that only 
a small range of frequencies need be avoided. Wall effects are shown 
theoretically to be more pronounced as the Mach number is increased and 
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at high Mach nvimhers are found to he large even at frequencies well 
removed from resonance. 


langley Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
langley Field, Va., January 12, 1955- 
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APPENDIX 

RTOUCTION OF IMTEGRAL EQUATION TO THE CASE OF ZERO FREQUENCY 


In this appendix, the integral equation for the downwash for a wing 
oscillating in a compressihle medium in the presence of wind-tunnel walls 
is reduced to the zero- frequency condition. 

If equation (l) of the text is written as 

w(x) = lim r L(xo)foiC(M,z) + oiC(M,z,H^dXo (A1) 

0)— >0 pu ^^-1 L . J 

and the limit taken as oi — > 0, it will he found that all the terms of 
oK(M,z) and oiC(M,z,H) vanish except terms involving Hq(2) , These 
terms become infinite; however, as o) — >0, the asymptotic expansion 
for very small values of the argument may be used. Therefore, 


Hi(2)(nRn) 


2 

itipRn 


and 

UM = W(x - Xq) 

o)->0 rt[(x - Xo)2 + 3^(nH)2] 

The vertical induced velocity may then be written as 


v(x) 



+ 2 ^ (-1)^ 


n=l 




X - Xq 


dxo 


(A2) 


or 


v(x) =- 


Nbb 

J.i 


— + 2 ^3 (-1)*^ 


Jt(X - Xq) 



■^X - Xq) n-1 

pH 


3%^ 


dxo 


(A5) 
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Equation (A5) may "be written as 



The additional induced velocity due to the presence of tunnel walls 
for the steady-state case in compressible flow is given hy equation (4o) 
of reference 15 . Equation (A2) can be reduced to the same form by making 
the approximation that the airfoil chord is small congjared with the tunnel 
hei^t. ■ 
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Figure 2 .- Schematic drawing of model and oscillating mechanism. 
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Figiire 3«- Schematic diagram of instrumentation. 
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Figure 5.- Effect of variation in Mach nvunber at constant height- 
semi chord ratio H = 7»60 on the magnitude and phase angle of 
the lift of a pitching wing. 
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Figiire 6.- Effect of variation in Mach n\mber at constant height- 
semi chord ratio H = 7*60 on the magnitude and phase angle of 
the lift of a translating wing. 
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Figiire 7-- Effect of variation in Mach number at constant height- 
semi chord ratio H = 7 *60 on the magnitude and phase angle of 
the moment of a pitching wing. 
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(b) Phase angle of the lift 














NACA TN 3416 


4l 



NACA- Langley - 6-13-55 - 1000 


Figure ^10.- Comparison of lift ratio for M = 0 and M = 0.3'. 



